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ABSTRACT 

We predict the performance of the Planck satelhte in determining the bulk flow through kinetic 
Sunyaev-Zeldovich (kSZ) measurements. As velocity tracers, we use ROSAT All-Sky Survey 
(RASS) clusters as well as expected cluster catalogs from the upcoming missions Planck and 
eRosita (All-Sky Survey: EASS). We implement a semi-analytical approach to simulate realistic 
Planck maps as well as Planck and eRosita cluster catalogs. We adopt an unbiased kinetic SZ 
filter (UF) and matched filter (MF) to maximize the cluster kSZ signal to noise ratio. We find 
that the use of Planck CMB maps in conjunction with the currently existing ROSAT cluster 
sample improves current upper limits on the bulk flow determination by a factor ^ 5 ('~ 10) 
when using the MF (UF). The accuracy of bulk flow measurement increases with the depth and 
abundance of the cluster sample: for an input bulk velocity of 500 km/s, the UF recovered velocity 
errors decrease from 94 km/s for RASS, to 73 km/s for Planck and to 24 km/s for EASS; while 
the systematic bias decreases from 44% for RASS, 5% for Planck, to 0% for EASS. The 95% 
upper limit for the recovered bulk flow direction Aa ranges between 4° and 60° depending on 
cluster sample and adopted filter. The kSZ dipole determination is mainly limited by the effects 
of thermal SZ (tSZ) emission in all cases but the one of EASS clusters analyzed with the unbiased 
filter. This fact makes the UF preferable to the MF when analyzing Planck maps. 

Subject headings: Cosmology: cosmic microwave background, observations, diffuse radiation 



1. Introduction 

Peculiar velocities, along with inhomogeneities, 
can be used to constrain cosmology. A coher- 
ent, large scale peculiar velocity, also called bulk 
flow, may originate from spatial inhomogeneities 
in the mass distribution of large scale structures 
around us. The standard inflationary model pre- 
dicts that the rms bulk velocity within a sphere 
of radius R decreases linearly with comoving dis- 
tance in the ACDM universe with Kms('' > 50 — 
100 Mpc h-i) « 250( ^°" '^P'^ '''' ) km/s (Kashlin- 
sky & Jones 1991). Galaxy cluster peculiar ve- 
locity surveys at scales i? < 60 Mpc h~^ gener- 
ally agree with theoretical predictions of the clus- 
ter bulk velocities. However, recent measurements 
at larger scales {R > 100 Mpc h~^) indicate that 
the bulk flow velocity is significantly larger than 
the ACDM prediction with statistical significance 



up to 3cr (Kashlinsky et al. 2008; Feldman et al. 
2010). 

At scales i? < 60 Mpc h~^, an enhancement 
of the bulk flow with respect to the predicted 
ACDM value in the local universe is attributed to 
a large scale void or overdensity at these depth. 
This is thought to be the cause of the Local 
Group (LG) motion with respect to the CMB 
rest frame, with velocity v — 627 ± 22 km/s to- 
wards I = 276°, b = -1-30°, in alignment with the 
CMB dipole (Kogut et al. 1993). The measured 
value is within the cosmic variance limit. Dressier 
et al. (1987) and Lynden-Bell et al. (1988) iden- 
tified the Great Attractor (GA), a mass concen- 
tration of ^ 10^^ Mq, as the origin of the fiow 
at < 60 Mpc h^^. On large scales where ACDM 
predicts bulk fiows with negligible amplitude, ob- 
servations show the contrary; Lauer & Postman 
(1994) found a strong bulk flow signature with 
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V = 561 ± 284 km/s towards I = 220°, b ^ -28° 
(±27°) at a depth of 110 Mpc h'^. Using a sam- 
ple of 119 Abell clusters within 150 Mpc h~^ they 
found that the flow originate from a mass con- 
centration beyond 100 Mpc h~^. Similarly, Feld- 
man et al. (2010) found a bulk flow on scales of 
^ 100 Mpc h-i with w = 416 ± 78 km/s towards 
I = 282° ±11°, 6 = +6° ±6°, in disagreement with 
the WMAP5 cosmological parameters at 99.5% 
confidence. The direction and scale of these bulk 
flows are shown in Figure 1. 

Several theories have been suggested to ex- 
plain the high bulk flow velocities at large scales. 
One explanation is that pre-inflationary fluctua- 
tions in scalar fields on superhorizon scales gives 
a titled universe (Turner 1991, Kashlinsky ct al. 
1994, Mersini-Houghton k Holman 2009). In this 
picture, matter slides from one side of our Hubble 
volume to the other, producing an intrinsic CMB 
dipole anisotropy as seen in the matter rest frame. 
This inhomogeneity generates a bulk flow with cor- 
relation length of order the horizon size. Alter- 
natively, Wyman & Khoury (2010) have showed 
that a strengthened gravitational attraction at late 
times can speed up structure formation and in- 
crease peculiar velocities. Using N-body simula- 
tions, they found an enhancement in large scale, 
R > 100 Mpc h"-'^, bulk flow velocities of up to 
~ 40% relative to the ACDM cosmology. A sim- 
ilar approach using modified gravity is discussed 
in Afshordi et al. (2009) and Khoury & Wyman 
(2009). 

Kashlinsky & Atrio-Barandela (2000) have pro- 
posed a method aimed at determining the largest 
scale bulk flows from galaxy cluster peculiar veloc- 
ities measured using the kinetic Sunyaev-Zeldovich 
(kSZ) effect. If many galaxy clusters are mov- 
ing with a coherent motion with respect to the 
CMB rest frame, the kinematic part of the SZ sig- 
nal acquires a dipole moment. Since the kSZ sig- 
nal is proportional to line of sight velocity, such 
a measurement directly probes the bulk flow, free 
of distance measurement errors. Several authors 
attempted to measure the bulk flow using the 
measured kSZ effect in WMAP data. Kashlinsky 
et al. (2008) (hereafter KAKE, and later Kashlin- 
sky et al. 2010) first utilized this method, claim- 
ing a large-scale flow with v > 600 km/s out to 
~ 575 Mpc h"-'^, without sign of convergence to 
the ACDM predicted value. However, by repeat- 
ing the same method Keisler (2009) did not de- 
tect a statistically significant bulk flow. Osborne 
et al. (2010) (hereafter OMCP), used filters con- 
structed to enhance the signal to noise of the ki- 



netic signal and found no significant velocity dipole 
in the WMAP 7 year data. More specifically, they 
found a 95% bulk flow upper limits of the order 
of 4600 km/s in the direction of the KAKE They 
also showed that the matched filter outperforms 
the unbiased one when WMAP data are used, 
and demonstrated that CMB and instrument noise 
dominate the uncertainties. 

In this work, we apply the scheme of OMCP 
to study the capability of Planck data and future 
cluster surveys to measure bulk flows. The use of 
Planck maps is expected to produce improved re- 
sults with respect to the WMAP case because of 
reduced instrument noise, wider frequency cover- 
age (ensuring better foregrounds' subtraction) and 
increased spatial resolution of this mission. In ad- 
dition, Planck will also produce the first all-sky SZ 
survey with a median redshift of z = 0.3 {Planck 
Blue Book) and will ensure a better performance 
in separating the tSZ from the kSZ signal for any 
cluster sample considered. 

We investigate to what extent the use of Planck 
maps, in combination with data for existing 
ROSAT clusters, improves on bulk flow determi- 
nation from WMAP. We also assess the expected 
performances of bulk flow measurements for up- 
coming all-sky cluster catalogs, such as the ones 
derived from Planck and cRosita satellites. Such 
samples are more abundant and extend to higher 
redshifts than the one in hand. 

Our goals are: (i) to determine the sensitivity 
of the cluster velocity dipole measurement with 
the Planck specifications and assess the nature of 
the uncertainty; (ii) to study which cluster survey 
can best constrain the bulk flow; (iii) to study the 
performance of the filters used in OMCP with the 
Planck setup. 

This paper is organized as follows. The bulk 
fiow velocity expected from the ACDM model is 
calculated in section 2. In section 3, we briefly de- 
scribe the procedure we use to measure the bulk 
flow velocity. In sections 4 and 5 we give details 
of the SZ and X-ray cluster catalogs we use and 
describe the procedure we adopt to generate sim- 
ulated SZ maps. In section 6, we present the two 
filters we use to reconstruct the kSZ signal from 
the CMB maps. In section 7, we describe the anal- 
ysis pipeline we use to measure and calibrate the 
cluster dipole. In section 8, we describe the sys- 
tematic effects that may contaminate our results. 
The results are presented in section 9, followed by 
our conclusions in section 10. Throughout this pa- 
per, we assume a ACDM cosmological model with 
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Fig. 1. — The dipole direction in Galactic coordinates of current bulk flow measurements, (a) CMB b = 
48.26 ± 0.03°, I = 263.99 ± 0.f4° (Jarosik et al. 2010); (b) Local Group b = 30±3°J = 276 ± 3° (Kogut 
et al. 1993); (c) KAKE b^34°, I ^ 267°; (d) Watkins et al. (2009) 6 = 8 ± 6°, Z = 287 ± 9°; (e) Lauer & 
Postman (1994) b = -28 ±27°, I = 220 ±27°; (f) Hudson et al. (2004) b = 0±11°, I = 263 ±13°; (f) Feldman 
et al. (2010) 6 6 ± 6°, Z = 282 ± 11°. The size of the colored region is proportional to the amplitude of the 
measured bulk flow. The color code represents the convergence depth of the bulk flow, in units of Mpc h~^. 
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flm = 0.3, flA = 0.7, h = 0.72, w = -1, ag = 0.8. 
2. Theoretical Aspects 

2.1. The Kinetic SZ effect as a probe of 
velocity 

The SZ effect (Sunyaev & Zeldovich 1970) is a 
secondary CMB anisotropy caused by the scatter- 
ing of CMB photons by high energy electrons, such 
as those inside the intra-cluster medium. The frac- 
tional change in temperature of the CMB photons 
is the sum of two componets: ST = ATksz+ATtsz, 
where the first term and second terms are the ki- 
netic (kSZ) and termal (tSZ) Sunyaev-Zeldovich 
effect respectively. The kSZ is caused by the 
Doppler shifting of CMB photons due to the pe- 
culiar motion of the galaxy cluster with respect to 
the CMB rest frame. The fractional temperature 
change is: 



\ TCMB / 



(1) 



kSZ 



where Vp is the linc-of-sight peculiar velocity of 
the cluster, t = J ricdl is the optical depth of 
the cluster, Ue is the electron density, and ctt is 
the Thompson scattering cross section. For typ- 
ical galaxy clusters with r ^ 10""^ and v ^ 500 
km/s, we expect a CMB temperature change of 
ATksz ^ 5fj,K at the location of a galaxy clus- 
ter. To date, the kSZ effect has not been mea- 
sured in individual galaxy cluster. Nevertheless, if 
many galaxy clusters are moving in the same direc- 
tion, their velocity field creates a dipolar pattern 
in the CMB radiation at large scale and leads to 



a net dipole moment Ci^ksz = ^cmb 



(KAKE), where (r) is the mean optical depth of 
the cluster sample. 

The thermal Sunyaev-Zeldovich (tSZ) effect is 
the boosting in energy of CMB photons scattered 
by hot electrons in the intra cluster medium and 
results in a fractional temperature change of: 



V-tCMB/, 



(2) 



tsz 



where y = (Tt J dl{k-B,Te) / {me(?)nf, is the Comp- 
tonization parameter, Tg is the electron tempera- 
ture of the cluster, f{x) = x(e^ -|- l)/(e'^ — 1) — 4 
is the frequency dependence of the tSZ effect and 
X = hu/kBTcMB- The minimum and maximum 
temperature changes occur at u = 143 GHz and 
u = 353 GHz with a null at = 217 GHz. For typ- 
ical clusters, Tg ~ 10 keV gives ATtsz ~ WOfiK. 



In individual clusters, the tSZ effect is typically 
larger than the kSZ effect: 
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Nevertheless, the dipole component of the kSZ ef- 
fect may dominate over the statistical dipole com- 
ponent of the tsz effect if a bulk flow is present. 
Furthermore, the different frequency dependence 
of the two effects allows us to extract the kSZ sig- 
nal. 

2.2. Expected Bulk Flow Velocity from 
ACDM model 

In the linear theory of structure formation the 
peculiar velocity field of galaxy clusters is related 
to the matter overdensity through the continuity 
equation: 



dPijz) 
dz 



(4) 



and fi is 



where f{z) ^ = H{z) 

the growth factor (equation B2). 

At scales much larger than the size of a cluster, 
the peculiar velocity field is correlated within a 
given region leading to a coherent motion, or bulk 
flow, of the objects inside the region. Assuming 
the distribution of clusters is isotropic, the bulk 
flow velocity within a region of size R is given by: 



dk 



Pmjk) 

27r2 



\W{kR)\' (5) 



where W{kR) is the Fourier transform of the top- 
hat window function, R = c/H{z')dz' is the 
comoving radius and Pm{k) is the present-day lin- 
ear matter power spectrum. The cosmology de- 
pendence is embedded in the matter power spec- 
trum Pm{k) and its time evolution. In particular 
the rms velocity, cr„, is proportional to the am- 
plitude of the matter power spectrum erg which 
is known with an uncertainty < 10% (e.g. Larson 
et al. 2010). 

We incorporate the effect of the different red- 
shift distributions of the clusters in the catalogs 
we use by weighting the bulk flow velocity by a se- 
lection function (t){z). The selection function takes 
into account the fact that the survey is not com- 
plete at all of the redshifts we use. Then, 
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where V{z) is the comoving volume, (f){z) is the co- 
moving number density or selection function. We 
calculate (j){z) from the halo mass function (de- 
scribed in appendix B.l) and cluster mass limits 
for each survey (section 4), 



(t){z) = n{z) = 

J A 



dM^A^ (7) 



where Mjnin{z) is the limiting mass of object in 
the cluster survey at redshift z. (j){z) is normalized 
such that 



/" 

Jo 



(t){z'){dV{z')/dz'dn)dz' = 1 



Figure 2 shows the rms bulk flow velocity pre- 
dicted by equation 6 for the three cluster samples 
considered in this work. The effect of the selec- 
tion function is small that the velocities among 
the three cluster samples are no different than 1% 
at all redshifts. For comparison, the rms bulk flow 
computed with no selection function is also plot- 
ted. The sample variance of the velocity distri- 
bution is the largest source of uncertainty. For 
a Gaussian density field the amplitude of the bulk 
flow has a Maxwellian distribution with the proba- 
bility of having a velocity lying between V and V+ 
dV given by P{V)dV oc V"^ expi^l.bV"^ /a^)dV . 
The 95% confldence limits on the measured bulk 
flow of amplitude V are then V/3 < V < 1.6V. 
This is the shaded region in the plot. 

In addition to large scale correlated cluster ve- 
locities, each cluster has a peculiar velocity caused 
by matter inhomogeniety on cluster scale (i? ~ 
8 Mpc h~^). We include this components in our 
simulations. 

3. Methodology 

3.1. Assumptions cind Definitions 

In this section we define the different physical 
processes contributing to the sky signal, as well as 
the Planck instrumentation. 

The observed dipole signal at clusters' locations 
has, in principle, several contributions: 
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where a\j-^ are the dipole coefficients of the CMB, 
instrument noise, thermal SZ, kinetic SZ, extra- 
galactic radio and infrared point sources, and 
galactic signal. An appropriate temperature map 
mask can be used to suppress the galactic com- 
ponents and flt for the cluster velocity dipole 
outside the galactic region. We therefore do not 
consider the effect of galactic foregrounds in this 
work and apply a cut at \h\ < 7° to the simu- 
lated Planck maps. The contribution from in- 
frared point sources is not well characterized at 
the present time and no model is available to ac- 
curately describe their abundances within galaxy 
clusters. We therefore do not consider the infrared 
point source signal in this work, and only simulate 
the first 4 terms in equation 8 as well as the radio 
point source term. 

Our goal is to separate a^f^ from the other 
components and calculate its amplitude in Planck 
maps using different cluster samples. The detector 
noise levels and beam sizes we use are summarized 
in Table 1. Although CMB observations by Planck 
cover 9 frequencies from 30 GHz to 857 GHz, we 
only consider the 6 channels between 44 GHz and 
353 GHz because of large potential foreground con- 
tamination outside this range. 

3.2. Outline of the Method 

The steps we take to (estimate the cluster kinetic 
SZ dipole velocity are: 

1. We simulate cluster catalogs using halo 
model and cluster self similar scaling rela- 
tions. We do so for three different cluster 
samples from the ROSAT All-Sky Survey, 
Planck, and eRosita All-Sky Survey. 

2. Full-sky maps of the simulated kinetic and 
thermal SZ are created in the 6 Planck fre- 
quency channels together with realizations of 
the CMB. These maps are convolved with 
Gaussian beams and detector noise is added 
using the properties given in Table 1. 

3. The maps are filtered to enhance the kinetic 
SZ signal while suppressing the CMB and in- 
strument noise terms and removing the ther- 
mal SZ signal. 

4. The filtered maps are used to calculate the 
dipole moment of the reconstructed kSZ sig- 
nal. The dipole is calculated at the cluster 
positions using different redshift shells. 

5. Temperature dipoles (in units of Kelvin) are 
converted into velocity dipoles (in units of 
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Fig. 2. — The theoretical bulk flow velocity in ACDM cosmology, using selection functions for the three 
cluster surveys, smoothed over a top hat window function W{kR) with the comoving sphere centered on 
the observer with radius R. Also plotted is the amplitude of the expected bulk flow with uniform selection 
function. The shaded region is the uncertainty of the expected bulk flow from sample variance. 



Table 1 

Characteristics of the Planck channels 



Planck Channel 


1 


2 


3 


4 


5 


6 


Center frequency vq (GHz) 
Resolution A0 (FWHM) 
Noise level ctn (77^, lO^^) 


44 
26. '8 
1.1 


70 
13.'1 
2.2 


100 
9. '2 
0.6 


143 
7.'1 
1.0 


217 
5'.0 
1.6 


353 
5.'0 
4.9 



Note. — Characteristics of the Planck iF/-receivers (column 2-3) 
and ffi^J-bolometers (column 4-7): center frequency z/q, angular res- 
olution AO in FWHM, and instrument noise variance per pixel ctn 
(thermodynamic temperature units). 
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km/s) by a calibration matrix M, which is 
calculated from kSZ realizations with pre- 
defined bulk flow amplitudes, i.e. ay = 
Max, where ax and ay arc the monopole 
and dipole coefficients in temperature and 
velocity units respectively. 

We iterate the above pipeline to study the prop- 
erties of the recovered bulk flows by (a) using dif- 
ferent filters; (b) varying the amplitude of the in- 
put bulk fiow velocity; (c) using different cluster 
catalogs; (d) inputing different systematic compo- 
nents into the maps to determine their importance. 

3.3. Optical depth 

The calibration of conversion from dipole am- 
plitude to flow velocity, through the calibration 
matrix M (section 7.2), requires cluster position 
and optical depth. In this work, we assume op- 
tical depth is measured in the observed samples 
with negligible measurement errors and consider 
the scatter in the Y — M relation for Planck and 
EASS clusters and L—T relation for RASS clusters 
when producing our simulations (appendix B.2). 
When dealing with clusters from observations, we 
have to derive the optical depth of the cluster sam- 
ple to reconstruct the kinetic SZ signal. We briefly 
outline below how we can recover optical depth 
from SZ and X-ray clusters respectively. 

For SZ observation, Planck, the optical depth 
of each cluster is directly obtained from their SZ 
flux, or equivalently the integrated compton Y- 
parameter, by t = Y {me(?) / {kBTe) if the elec- 
tron temperature Tg is known. The tempera- 
ture can be obtained from X-ray measurements 
through the L — T scaling relations. For X-ray 
observations, RASS and EASS, we can calculate 
the optical depth from the electron density by 
r = (Tt / n^dl. We follow OMCP to determine 
the cluster electron density (see equation 8 in 
OMCP) using the Bremsstrahlung emission lumi- 
nosity (Gronenschild & Mewe 1978) and the re- 
lation Tx = (2.76 ± 0.08)L^^^=^°-°^ (White et al. 
1997) to calculate the electron temperature of the 
cluster. For the RASS clusters, propagating the 
error in observed Lx to optical depth gives uncer- 
tainty of - 15%. The error in Lx for the EASS 
clusters is likely to be smaller since the backgroimd 
noise for eRosita is expected to be lower. The un- 
certainty in the optical depth for RASS is therefore 
an upper limit for EASS clusters. 



4. Cluster Catalogs 

We construct cluster catalogs for three represen- 
tative all sky surveys: the combined ROSAT All- 
Sky Survey (hereafter RASS), the future Planck 
cluster survey and eRosita All-Sky Survey (here- 
after EASS)^. Each survey has its own advantages: 
the RASS sample already exists and has been used 
for previous bulk flow measurements (e.g. KAKE 
and OMCP ); the Planck cluster sample will con- 
tain about three-flve times as many clusters, will 
extend to greater distances and, being an SZ sur- 
vey, will be subject to different selection effects 
than RASS; the EASS cluster sample will be larger 
than both the ROSAT and Planck samples and 
will extend to greater redshift z ~ 1. All three 
surveys provide isotropic samples outside regions 
at low galactic latitudes. All surveys are fiux- 
limited, and so at any redshift z only objects of 
mass M > Miim(-Piim) are included in the catalog. 
A summary of the catalog properties is listed in 
Table 2. 

All surveys are flux-limited, and so at any red- 
shift z only objects of mass M > Miini(Fiiiji) are 
included in the catalog. A summary of the catalog 
properties is listed in Table 2. 

4.1. eRosita All-Sky Survey 

The eRosita mission^ is expected to be launched 
in 2012 and perform the first all-sky X-ray imag- 
ing survey in the X-ray energy range up to 5 
keV with a limiting flux of i^um.o.s-s keV = 1-6 x 
10-13 erg s-i cm-2. The EASS is expected to 
yield a catalog of a few tens of thousands of clus- 
ters out to redshift z « 1.3. Multi-band optical 
surveys are planned to provide photometric and 
spectroscopic redshifts for the EASS clusters with 
masses above 3.5 x IO^^'Mq/i-^ (Cappelluti et al. 
2010). The second phase of the eRosita survey 
(the wide survey) will detect more clusters by us- 
ing longer exposures than the flrst phase (the all 
sky survey). However, it will only cover about half 
of the sky resulting in a non-isotropic cluster sam- 
ple and so we only use clusters expected from all 
sky survey. We apply a galactic cut at \b\ < 7° 
giving a sky fraction 0.727. The limiting mass of 
cluster to be included in the sample at redshift z 
is given by (FedeH et al. 2009): 



^ http;//www.mpe.mpg.de/projects.html#erosita 
^mission definition document at 

http://www.mpe.mpg.de/erosita/MDD-6.pdf 
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Fig. 3. — The minimum mass as a fmiction of redshift for a cluster to be included in the RASS, Planck, 
and EASS catalogs. The cluster mass that we use in the catalogs is max[lO^'*M0, Miii„(2)]. For the RASS 
sample, we fit equation 9 to the archived catalog and obtain the values shown in this plot. 




Fig. 4. — The differential cluster redshift distribution per square degree for clusters in the Planck, EASS 
and RASS, with cluster mass M > max[lO^'*M0, Mii„i(2)]. The ACDM cosmology is assumed and as — 0.8 
is used for the number counts of the Planck and EASS clusters. The observed number count of the RASS 
clusters are overlaid with the solid lines for comparison with the fitted number counts (dash dot line) . 
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im,200 
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E{z) 



lim.0.5-2kcV 



/Cb 



1/1 



1.097 X lO^Sergs^i 
(9) 

where E{z) = H{z)/Hq is the normalized Hub- 
ble parameter, is the luminosity distance 
and Cb is the band correction factor which converts 
the bolometric flux to the eRosita energy 0.5 — 2 
keV. We estimate Cft by assuming a Raymond- 
Smith (Raymond & Smith 1977) plasma model 
with metalliticity of OAZq, a cluster temperature 
of 4 keV, and a Galactic absorption column den- 
sity of riH = 10^^ cm~^. For consistency, we use 
the virial mass definition throughout this work and 
we convert M200 to My (and write My as M here- 
after) using the conversion fitting formula by Hu 
& Kravtsov (2003). 

To create mock catalogs, we assign a mass and 
redshift to each cluster in the EASS sample and 
simulate the properties of the clusters using X-ray 
scaling relations (described in appendix B.2). 

4.2. The Planck Catalog 

Planck is imaging the whole sky with an un- 
precedented combination of sensitivity (AT/T ~ 
1 X 10~^ at 143 GHz), angular resolution (5' at 
217 GHz), and frequency coverage (30 — 857 GHz). 
The SZ signal is expected to be measured from 
a few thousand galaxy clusters. Planck will pro- 
duce a cluster sample with median redshift ~ 0.3 
with significant fraction beyond redshift 1 (Planck 
Blue Book). As with EASS, we do not simu- 
late clusters with galactic latitude \h\ < 7° to 
minimize the Galactic signal. The SZ observ- 
able is the integrated Comptonization parameter 
Y = J y dficiustcr (appendix B.2) with expected 
values of F200 > 10~^arcmin^ (Malte Schafer & 
Bartelmann 2007), where I200 is the integrated 
comptonization parameter within r2oo • Fedeli et al. 
(2009) provide a fitting formula for the limiting 
mass as a function of redshift based on simulations 
by Malte Schafer & Bartelmann (2007): 



lim,200 



(^) 



1015Mq/i 



g-1.200+1.469tan-i[(2-0.11)]'' jf ^ > 
g-1.924+8.333. jf ^ < 0.11 



(10) 



4.3. ROSAT All-Sky cluster survey 



The RASS sample, consisting of clusters from 
the REFLEX (Bohringer et al. 2004), BCS (Ebel- 
ing et al. 2000b), eBCS (Ebeling et al. 2000a), 
CIZA (Ebeling et al. 2002; Kocevski et al. 2007), 



and MACS (Ebchng et al. 2001) catalogs com- 
554 prises a total of 827 clusters. All of the clusters 
have spectroscopic redshifts with z < 0.5. We sim- 
ulate the SZ signal of the clusters in this catalog 
from the observed X-ray properties via scaling re- 
lations, while the redshifts and positions on the 
sky are taken from the catalog. 

Figure 3 shows the limiting mass for clusters to 
be included in the RASS catalog. We fit equa- 
tion 9 to the 827 RASS clusters and find an effec- 
tive fiux limit of Jiim = 1.2 x 10~^^ergs~^ cm~^. 
In Figure 4 we show the cluster number counts 
dN/dz/dO,, for our cluster catalogs. 

5. Sky Simulations 

We create full-sky SZ maps using a semi- 
analytical approach. We follow the method pre- 
sented in in Delabrouille et al. (2002) and Waiz- 
mann & Bartelmann (2009) that used the halo 
model and cluster gas properties. The simulations 
include SZ emission, primary CMB anisotropy 
and instrument noise, while diffuse Galactic fore- 
grounds are excluded. 

Clusters arc simulated with properties expected 
from the surveys described in section 4 with num- 
ber densities given by the halo mass function 
N{z,M). Each simulated cluster is then assigned 
a mass and redshift and the gas properties are then 
derived. We leave the details of the number counts 
and the cluster models to appendix B and focus on 
the preparation of the SZ maps here. 

In the nonrelativistic limit the distortion to the 
CMB temperature of the kinetic SZ effect is given 
by equation 1 but with the peculiar velocity re- 
placed by iT = T4,uik + ^^pecuiiar- We give an over- 
all bulk velocity Vhu\k to the whole cluster sam- 
ple which is a free parameter of the simulation. 
Besides a bulk motion, we also give each cluster 
a random peculiar velocity ^peculiar drawn from a 
Gaussian distribution with variance given by equa- 
tion 5 with R = 8 Mpc h~^ 

Maps of CMB anisotropics are generated from 
the angular power spectrum using the CMBfast 
P ^^<|:ode (Seljak & Zaldarriaga 1996) for a fiat ACDM 
iosmology. The CMB and SZ maps are com- 
bined and then smoothed with a Gaussian beam 
with sizes given in Table 1. Finally, the instru- 
ment noise is added with the noise variance and 
the beam size of the Planck channels are listed 
in Table 1. Synthesized maps are generated in the 
HEALPix pixelization scheme (Gorski et al. 2005), 
with resolution of Nside=1024. 
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Table 2 

Parameters of cluster catalogs 



Catalog 


RASS 


Planck 


EASS 


units 


/o.5-2keV 


Y200 


/o.5-2keV 




ergs^^cm~'^ 


arcmin^ 


ergs~^cm~'^ 


Flux limit 


1.2 X 10-1=^ 


10"^ 


1.6 X 10"^^ 


z* 

median 


0.091 


0.246 


0.211 




827 


2700 


33000 


(r) 


(7.1 ±2.2) X 10-3 


(8.5 ± 1.7) X 10-3 


(6.1 ± 1.5) X 10-3 



Note. — * The median redshift is computed from the simulated cluster 
sample with a cut-off redshift z=0.5 for all three cluster samples. 



6. Filtering: Reconstruction of the Kinetic 

SZ signal 

The kinetic SZ signal from a cluster is embedded 
in the CMB and instrument noise, and tSZ emis- 
sion is also present at the same location. We there- 
fore filter the maps to increase the signal-to-noise 
of our cluster dipole measurement. We use two 
types of linear multifrequency filters to reconstruct 
the cluster kinetic SZ signal. The filters are con- 
structed with the aim of minimizing the CMB and 
noise variance in the map. We leave the derivation 
of the filter shapes to the appendix A. The first fil- 
ter is a matched filter (hereafter MF; Haehnelt & 
Tegmark 1996; Herranz et al. 2002) that is opti- 
mized to detect the kSZ signal. The second filter 
(hereafter UF) is subject to the additional con- 
straint of removing the tSZ signal at the cluster 
location. The UF was first proposed by Herranz 
et al. (2005) for use on flat patches of the sky and 
later Schafer et al. (2006) extended the scheme for 
use on full sky maps. Both filters use all of the fre- 
quency channels and take into account the statis- 
tical correlation of the CMB and instrument noise 
between different frequencies. 

The filter kernels are: 

c 

matched filter : = ^— Bi (11) 

7 

unbiased filter : = ^(aBi - ^Fi) (12) 

where a = JZlZo FTCf^F,, /3 = E!=r B^'Cf^Fi, 
7 = E!=r Bi^Cr^Bi, and A = aj- p\ 

We assume the clusters are point sources in 
the Planck maps and so Bi are the spherical har- 



monic coefficients of the Planck beam function, 
Fi = /(i^)Bi where f{u) gives the tSZ frequency 
dependence, Ci = Cp'^^ -|- Cf°'«^ is a matrix 
at every multipole giving the sum of the cross 
power spectra of the CMB and instrument noise 
between frequency channels. The filters are nor- 
malized such that the filtered field gives the am- 
plitude of the kinetic SZ signal at the central pixel 
of each cluster. Given the spatial resolution of 
Planck we compute the filter kernels out to multi- 
pole /max = 3000. Figure 5 shows the MF and UF 
for the 6 Planck frequency channels we use, using 
the Planck beam FWHM and pixel noise variance. 

6.1. Matched filter (MF) 

The MF is less complicated than the UF since 
most of the features are only seen at the cluster 
scales (1000 < I < 2500) while the smallest and 
largest scales are suppressed. As with the UF, the 
channels v = 100, v = 143, p = 217 GHz are given 
the largest weights. The matched filter is more 
efficient than the UF at removing the CMB and 
instrument noise components because it forces the 
scales at which the noise and CMB dominant to 
disappear. However, since the MF is not designed 
to remove the thermal SZ signal, the amplitudes 
of the filter kernels are all positive. 

6.2. Unbiased Kinetic SZ filter (UF) 

The UF is constructed to give an unbiased esti- 
mate of the kinetic SZ signal at the clusters' loca- 
tions. At large and intermediate scales (l < 1000) 
the channels below 1/ = 100 GHz are subtracted 
from those above. The CMB fluctuations which 
dominate the signal at these scales are therefore 
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1000 



2000 
multipoles I 



3000 



suppressed. At cluster scales (1000 < / < 2500) 
all of the filter kernels have positive amplitudes 
and the 100, 143, 217 and 353 GHz channels are 
given the largest weight. At 100 GHz and 143 
GHz the tSZ signals are negative, approximately 
zero at 217 GHz and positive at 353 GHz. The 
channels are weighted so that the tSZ signal is 
zero. Although the tSZ signal is smaller at 100 
GHz than the other three channels, the pixel noise 
level is lower (ctn ~ Q^K) and so the 100 GHz 
channel also has substantial weight. At smaller 
scales {I > 2500) the signal is dominated by in- 
strument noise and so all channels are suppressed 
by the filter. 

7. Analysis Pipeline 

7.1. Weighted Least Square Fitting of 
Dipole 

We fit the real spherical harmonic coefficients of 
4000 monopole, oq, and dipole terms, aix^ aiy and 
ai2, to the ffitered maps using a weighted least 
square fit which is based on the Healpix IDL pro- 
cedure remove_dipole (Gorski et al. 2005): 




ax = (X'^WX)-^X'^Wu 



(13) 



1000 



2000 
multipoles 



3000 



Fig. 5. — (Upper) Multi- frequency matched filter 
signed to detect the kSZ signal. (Lower) Unbi- 
ased multi-frequency matched filters designed to 
remove the tSZ components . Both filters are 
constructed using the Planck beam profiles as the 
cluster radial profiles. 



where ax is a vector of best fit monopole and 
dipole coefficients, u is the filtered map (e.g. equation A2), 
W is a matrix with diagonal elements equal to the 
weight given to each pixel of the map, and X is a 
matrix giving the contribution of the fitting func- 
tion to each pixel. We give the central cluster pixel 
of the map a weight Wi = 1 /cr^ j, where aN,i is the 
i~th pixel noise variance, all other pixels are given 
zero weight. The noise variances are calculated 
4000 from 100 filtered CMB and noise realizations. We 
also tried weighing each pixel by an estimate of 
the signal to noise, i.e. Wi = n/aj^ -, but we find 
a larger systematic bias due to tSZ contamina- 
tion in the recovered bulk fiow than when using 
the former weighting scheme. We therefore use 
a weighting scheme that involves only the pixel 
noise. 

We fit the dipole only at central cluster pixels 
because our filters are optimized to reconstruct the 
source amplitude if the source were centered at 
that pixel. 

7.2. Translating from fiK to km/s 

We construct a 4 x 4 matrix ay = Max to con- 
vert the dipole from temperature units (fJ-K) to ve- 
locity units (km/s), in which the matrix elements 
having units of km/s/K. This involves the use of 
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simulated kinetic SZ maps that ean aecount for 
the attenuation of the optical depth by the beam 
convolution and the filtering process. In principle 
one can perform this unit conversion in observa- 
tions, e.g. Planck temperature sky maps, using 
calibrated SZ simulations based on estimation of 
the optical depth of each clusters. This is the case 
in OMCP that they used the simulated kSZ maps 
of the RASS clusters to calibrate WMAP data. 
On the other hand, in this work we assume the 
optical depth, hence the kSZ signal, of each clus- 
ter is known and calculate the matrix M from our 
sky simulations. Note that the matrix elements de- 
pend implicitly on the average optical depth of the 
cluster sample, therefore it is specific to individual 
cluster experiment and we construct it separately 
for the three cluster surveys. We now describe how 
we create this matrix. 

We generate four different sets of kinetic SZ sig- 
nal only realizations with a given bulk flow ampli- 
tude: one with a monopole velocity and the other 
three with a dipole velocities in the x, y, z direc- 
tions. We choose an amplitude of of 10,000 km/s, 
which is large enough such that the recovered sig- 
nal is not confused by uncertainty of optical depth 
and random peculiar velocity. Since the kinetic SZ 
signal of each cluster is proportional to the optical 
depth, the elements of the matrix M depend im- 
plicitly on the average optical depth of the cluster 
sample. Each set of maps are passed through the 
pipeline to obtain the aim- The four aim from each 
set of realizations are the elements of each row of 
M~^. By repeating this procedure for the four ki- 
netic SZ signal only realizations, we construct the 
16 elements of the calibration matrix. We check 
that the off-diagonal elements of M are at most 
1% of the diagonal ones for all of the cluster sam- 
ples considered. 

We perform this exercise 20 times and assign 
the average value from these 20 realizations to the 
matrix elements. For each of the RASS, Planck, 
and EASS sample, we calculate one average M and 
apply it to the filtered temperature maps. Due to 
the intrinsic scatter of the Y — M relation intro- 
duced in the simulations, the average optical depth 
of any simulated cluster sample might deviate from 
the expected value from the characteristics associ- 
ated with that catalog, this may lead to calibration 
error (section 8.3). 

8. Error estimation 

The error budget to the dipole coefficients con- 
tains several terms: 



2 2 2 2 2 / \ 

f^av = «^CMB-Fnoise + "^kSZ + <^tSZ + <^PS (l^) 

where crcMB-i-noisc is the residual from the CMB 
and instrument noise, cr^gz is the error in the ki- 
netic SZ signal due to the intrinsic scatter of the 
Y — M relation and the random component of the 
galaxy cluster peculiar velocity that is not part of 
the bulk flow, a^g^ is the thermal SZ residual, and 
(Tpg is the contamination from radio or infrared 
point sources associated with clusters. In princi- 
ple, the last three may be correlated, while we are 
assuming here that correlations are negligible. 

The errors on the three dipole coefficients aim 
are correlated and the uncertainty in the best fit 
value can be described by a covariance matrix 
N. We compute the covariance matrix of each 
type of error in equation 14 by passing simula- 
tions of the components through our pipeline and 
performing the dipole fit on them. The scatter 
cr in dipole coefficients then provides an estimate 
of the noise correlations between the dipole direc- 
tions, i.e. N = (ayay). Then, 

= (av.rec - av,in)'^N~t (av,rec - ay.in) (15) 

where ay.rec is the best fit dipole coefficients of the 
recovered velocity and ay, in is the input velocity 
dipole coefficients. The confidence level of a given 
ay can be calculated from the probability dis- 
tribution with 3 degrees of freedom. 

Figure 6 shows the contribution from each sys- 
tematic effect to the dipole velocity. We discuss 
these results in the following subsections. 

8.1. CMB and instrument noise 

The CMB signal is not completely removed by 
our filters because there is CMB emission on clus- 
ter scales with the same frequency dependence as 
the kSZ signal. We evaluate the error from the 
CMB and instrument noise by measuring the bulk 
flow of 300 filtered CMB plus noise maps. As a 
further check on this procedure we also calculate 
the dipole on a single CMB plus noise realization 
and repeat the procedure with 100 realizations of 
the cluster sample, each having a different spatial 
distribution. The results from the two procedures 
are similar with errors that are both gaussian with 
similar variance. 

As expected the UF errors are larger than the 
MF ones by about 30% for the three cluster sam- 
ples. We find that the error is largest for the RASS 
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RASS 




1 1 1 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 

< Z < Z 

Planck 




1 1 1 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 

< z < z 

EASS 




1 1 1 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 
< Z < Z 



Fig. 6. — The systematic error contributions to 
the bulk flow measurement using the RASS (top), 
Planck (middle), and EASS (bottom) cluster sam- 
ples: instrument noise and CMB (cross), thermal 
SZ (asterisk), uncertainty in the kinetic SZ sig- 
nal due to scatter in both the optical depth and 
cluster peculiar velocity (diamond), and the con- 
tamination from radio point sources (triangle for 
our simulation and solid line for PSM simulation). 
The data points are the average values from 100 ^3 
realizations and the error bars are the standard 
deviations of the one-sided distribution. 



sample and smallest for the EASS sample, since 
the EASS sample contains more clusters. 

8.2. Thermal SZ signal 

We expect the thermal SZ signal to be entirely 
removed when using the unbiased matched filter 
(UF). However, maps filtered with the MF suffer 
from a systematic bias by the unfiltered thermal 
SZ signal that contaminates the kinetic SZ signal. 
In addition to the effect of filtering, the intrinsic 
scatter of the cluster y-parameter would also in- 
troduce further incomplete removal of the thermal 
signal and resulting in scatters around the mean 
recovered velocities. To evaluate the level of ther- 
mal SZ bias and uncertainty, we generate simu- 
lated maps containing only tSZ signal. 

We find that in thermal SZ signal only maps 
the average systematic bias (ilcr uncertainty as 
represented by the error bars in Figure 6) to the 
cluster dipole are v — 313 ± 48 km/s for RASS, 
V ~ 77 ± 33 km/s for Planck, and w 19 ± 8 km/s 
for EASS when filtered by the UF, and v = 885±66 
km/s for RASS, v ^ 171 ± 77 km/s for Planck, 
and V = 47 ± 20 km/s for EASS when fihered by 
the MF. The thermal bias is most serious for MF 
filtered RASS clusters while the la uncertainty is 
the largest for MF filtered Planck clusters. This 
is due to the fact that the average optical depth 
is the largest for Planck clusters and hence the 
largest intrinsic scatter (and uncertainty). 

We find a significant monopole: Vag = 730 
(RASS), 1440 (Planck), and 870 (EASS) km/s 
when fihered by the UF, and Va^ = 3500 (RASS), 
3890 (Planck), 2180 (EASS) km/s when filtered 
by the MF. The UF is more effective at removing 
the thermal SZ signal than the MF, by a factor 
of ~ 5. Though the monopole velocities are large 
compared to the thermal dipole, they have no ef- 
fect on the bulk flow measurement. 

8.3. Uncertainty in optical depth and ran- 
dom cluster peculiar velocity 

We find that the cluster random peculiar ve- 
locity not part of the bulk flow gives a negligible 
dipole with an average magnitude of w < 50 km/s 
for all three cluster samples and is independent of 
the type of filter used. The error due to uncer- 
tainty in the optical depth depends on the ampli- 
tude of the peculiar velocity since Cksz var ■ We 
find that ar ~ 15% for the Planck cluster sam- 
ple and ar < 15% (OlO) for the RASS and EASS 
clusters. The uncertainty for Planck clusters is 
relatively larger because it contains more abun- 



dant massive clusters and hence propagate larger 
errors to the optical depth for the same scatter 
in the cluster Y-parameter. As mentioned ear- 
lier, this uncertainty introduces a calibration er- 
ror when translating from /xK to km/s because of 
the difference of the average optical depth and the 
calibration matrix M. However, since we choose a 
large calibration velocity of 10,000 km/s, the cali- 
bration error is only about 5% and negligible when 
compared to the error due to CMB and noise. 

8.4. Extragalactic Point Sources 

At small scales extragalactic point sources are 

significant contaminants of CMB maps. They give 
a Poisson noise contribution to the measured angu- 
lar power spectrum and a non-Gaussian signature 
in the maps (Colombo & Pierpaoli 2010). Radio 
point sources have a falling spectrum with a ~ 0.7 
where S oc y~", whereas infrared point sources 
have a rising spectrum with typical spectral in- 
dices a ~ 3. Thus infrared sources tend to be 
the more significant contaminants at higher fre- 
quencies where the SZ effect is large. Recent mea- 
surements such as the South Pole Telescope (SPT) 
detected 47 IR sources at S/N> 4.5 (Vieira et al. 
2010). Nevertheless, the infrared point source sig- 
nal is less well characterized than the radio point 
source signal and reliable models are not available 
to estimate the number counts within a galaxy 
cluster. Furthermore, the 353 GHz channel, which 
suffers the most from IR point source emission is 
given less weight than the other channels at clus- 
ter scales. We therefore do not consider the contri- 
bution from infrared sources in this work. Radio 
sources are an important contaminant of SZ mea- 
siirements as they fill in the SZ temperature decre- 
ment if the soiirces are within the cluster being 
observed, leading to an underestimation of the SZ 
signal. We expect a negligible contribution since 
the cluster signal is strongest at frequencies where 
the radio point source signal is small. 

To estimate the level of radio contamination to 
the cluster kSZ dipole we simulate maps of the 
radio point source emission. To do this we use 
the radio luminosity function at 1.4 GHz of Lin 
& Mohr (2007) extrapolated to higher frequencies 
using, 

dlogP, J dlogFi.4 ■'^ ' ^ ' 

where Pi, is the radio luminosity in units of 
WHz~^at frequency v, a is the spectral index, 



and f{a) is the probability distribution of the 
spectral indices at 1.4 GHz, which is taken to be 
a Gaussian distribution with mean a = 0.51 and 
rms (Tq = 0.54. /(a) is taken to be f(a + 0.5) 
for i/ > 90 GHz to account for possible steepening 
of a at high frequencies. The RLF at 1.4 GHz is 
given by the fitting formula: 

(17) 

where u = 37.97, b = 2.40, x = 25.80, and 
u; = 0.78 (Lin & Mohr 2007). Here n is the radio 
source number density within r2oo- We integrate 
equation 16 to obtain the expected number of ra- 
dio sources, iVradio, in a cluster and assign fluxes 
to the A'radio sources using the RLF. Since bright 
sources will be masked in the Planck maps we also 
mask sources with S > where S'J,™ is the 

upper limit of the radio flux at frequency given 
by Lopez-Caniego et al. (2006, Table 2) leaving a 
radio map with sources having 5* < 500m Jy. 

We find no significant cluster dipole from radio 
point sources in the RASS, Planck, or EASS clus- 
ter samples. In the redshift shell z = 0—0.5 we find 
signals of < 10 km/s for EASS and v ^ 20 km/s 
for both Planck and RASS clusters when filtered 
with cither the UF or MF. We can further verify 
the result for the RASS sample by using more real- 
istic radio point source simulations that considers 
information from sources observed by NVSS, and 
simulations from the Planck Sky Model^ (PSM). 
As before, we mask bright radio sources in the 
PSM simulation. The result from the PSM sim- 
ulation is shown in Figure 6. The radio dipoles 
from the PSM are consistently larger than our ra- 
dio simulations at all redshift shells, but are of a 
similar order of magnitude with w ~ 50 km/s. 

9. Results 

We generate sets of simulations containing 
CMB, noise, thermal SZ signal, and kinetic SZ 
signal with 68 bulk flow velocities logarithmically 
spaced between 100 and 10,000 km/s. In this sec- 
tion we describe how well we recover the input 
bulk flow. 



^http : //www . ape . luiiv-par is7 . f r/APC.CS/Recherche/ 
Adamis/PSM/ psky- en . php 
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Fig. 7. — Recovered bulk flow velocity for clusters z < 0.5 ( Top: RASS, Middle: Planck, Bottom: EASS) 
from simulated maps with input bulk flow velocities logarithmically spaced between 100 and 10,000 km/s, 
filtered with the MF (black cross) and UF (blue diamond). The solid line represents perfect recovery. The 
dash lines quote the 95% confidence interval of the input bulk flow values. All results are based on maps 
containing kinetic and thermal SZ signals, but simulated with different systematic components. Plots labelled 
Noise and CMB are from maps containing instrument noise and CMB; plots labelled thermal SZ are from 
results containing thermal SZ emission and hence estimate the level of thermal SZ bias. 
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Fig. 8. — Recovered bulk flow velocity for three input bulk flows v =500 (black), 1000 (yellow), and 5000 
(cyan) km/s as a function of redshift for the Planck cluster sample. The recovered bulk flows are measured 
from simulations that contain all systematic errors we consider. 



9.1. Precision and detection threshold of 
Bulk Flow Measurement 

Figure 7 shows the recovered bulk flow values 
Kocovcrod versus the input values Mnput for all clus- 
ters with redshift z < 0.5. The dashed lines in 
the figures are the 95% upper limits to the re- 
covered bulk flow velocities and are computed by 
generating 10,000 Gaussian distributed dipole co- 
efficients at each of the x, y, and z directions 
having mean equal to the input bulk flow veloc- 
ity and thermal bias, and variance equal to the 
noise levels described in section 8, i.e. af™ = 
(ai"^ + a^l^"") ± a5'°f'=. The resulting magnitude 
of the dipole velocity follows the distribution, 
the magnitude and angle error (section 9.2) at 95% 
are then computed. 

We determine the precision of the velocity mea- 
surement by computing the deviation of the val- 
ues Vrecovered from their corresponding T^nput with 
a parameter defined to quantify the deviations as 
follows: 

/ Silo(log K-ecovered.i " log Vi„put,i)^ 

'^los^ = V 

(18) 

where the summation runs over the range of T^nput 
being considered and iV = 68. This parameter al- 
lows us to compare the level of dispersion among 
the results from the various experiments and sur- 



veys. We provide the values for each cluster sample 
in Table 3. 

Taking a representative velocity of l^nput =500 
km/s we find that the recovered velocities are over- 
estimated by 120% with a scatter of 16% when 
evaluated using the RASS clusters using the MF 
(^recovered = 1085 ± 82 km/s), but perfectly re- 
covered with only a scatter of 5% when evaluated 
on the EASS clusters using the UF (l^ccovcrcd — 
500 ± 24). The large overestimation in MF filtered 
RASS clusters is due to tSZ contaminations which 
has the largest effect for the RASS clusters. We 
find that the recovered velocity has the largest un- 
certainty for MF filtered Planck clusters, with a 
scatter of 25%, and is consistent with the error 
analysis that the uncertainty due to intrinsic scat- 
ter in Y-parameter is largest for Planck clusters. 
With the sensitivities we find, measurements using 
either the Planck or EASS cluster sample will be 
able to detect a bulk flow as large as that claimed 
by KAKE. If the bulk flow is consistent with the 
ACDM prediction of w = 30 km/s for redshift shell 
extending to z = 0.5, our analysis pipeline would 
recover a 95% upper limit to the velocity of, when 
the UF is used, v = 470 km/s for RASS cluster, 
f = 160 km/s for Planck cluster sample and v = 60 
km/s for EASS cluster sample. 

Figure 8 shows the recovered signal in different 
redshift shells at three input bulk flow velocities: 
Mnput = 500, 1000, and 5000 km/s, estimated from 
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100 realizations. For illustration purpose, only the 
result from the Planck cluster sample is shown. At 
these velocities, the recovered velocities are shown 
to be consistent with the input values within Icr 
at redshifts z = 0.1. At lower redshifts, the recov- 
ered velocities are noise dominated. At this red- 
shift limit, the Planck cluster catalog will contain 
approximately 400 clusters. If a bulk flow of am- 
plitude w 500 km/s is present our results suggest 
that a sample of about 400 clusters would allow 
Planck observations to constrain the bulk flow. 

9.2. Recovering the Bulk Flow Direction 

The direction of any bulk flow is not predicted 
by the ACDM model. We can estimate the er- 
ror in the direction of a detected dipole for each 
of the three surveys. For T^nput = 500 km/s to- 
wards I = 280° and b = 30° the average values of 
the dipole components, ((14), {Vy), {Vz)) from 100 
simulations, are listed in Table 3. The direction 
is recovered with deviation < 20% for the Planck 
and EASS cluster samples (except for (14) of the 
Planck sample filtered with the MF), but is un- 
constrained for RASS clusters. We find that the 
error is larger in the T4 and Vy directions due to 
the galactic cut. We also calculate the angle errors 
for a range of input velocities for four different in- 
put directions: the three perpendicular directions 
X (; = 0°, 6 = 0°), y (; = 90°, b = 0°) and z 
{I = 0°, b = 90°), as well as the reference Vlnput 
used above. The 95% upper limit to the angle er- 
rors. Attgs, is shown in figure 9. 

For maps filtered with the MF (UF) the 95% 
confidence limit on the dipole direction for a 500 
km/s input bulk flow is Aa|g° = 62°(34°) for 
RASS clusters, 34° (14°) for Planck clusters, and 
9° (4°) for EASS clusters. The errors are smaller 
than the discrepancies in the observed hulk flow di- 
rections. Therefore Planck should be able to bet- 
ter constrain the region where the bulk flow points 
to. 

9.3. Effect of the systematic error compo- 
nents on the amplitude of the recov- 
ered bulk flow 

The dominant systematic error is tSZ emis- 
sion, and the dominant statistical error is CMB 
and instrument noise and uncertainty in the Y- 
parameter. They are velocity-independent and 
dominates for v < 1000 km/s. At large velocities, 
the calibration error due to the scatter in optical 
depth dominates and gives a < 10% error in the 
velocity. 



Using inverse variance weights in the dipole fit- 
ting (as opposed to weighting by optical depth) is 
found to minimize the tSZ contribution. The re- 
covered bulk flows are generally overestimated by 
« 20% when W = r/a^^^^^ is used to weight pixels 
in the dipole fit. 

9.4. Performance of Filters 

The overall performance of both the MF and 
UF agree with expectations that maps filtered by 
the MF have higher kSZ signal to noise ratios while 
maps filtered by the UF have smaller thermal SZ 
bias. The UF is slightly a better filter in the Planck 
data since the dominant uncertainty of the mea- 
surement comes from the thermal SZ bias, and 
the recovered velocities are more accurately deter- 
mined. Nevertheless, the thermal SZ bias is not 
completely removed by the UF. This is because 
we assume a /3 model for the cluster profile, but 
assume a point source model to construct the fil- 
ters. We could assume a (3 model to construct the 
filters but we do not consider this here. 

We now compare the filter performance in the 
WMAP channels which are studied in OMCP (see 
their Figure 11) with the Planck channels in this 
work (Figure 7) , with the use of the RASS cluster 
sample which is used in both analysis. In OMCP 
the MF was found to be more sensitive than the 
UF, but we find the opposite. The velocity de- 
tected at 95% CL with the MF (UF) is reduced 
from ~ 5000 km/s (~ 10, 000 km/s) in the WMAP 
channels to ^ 1000 km/s (^ 1000 km/s) in the 
Planck channels. Considering that the amplitude 
of observed bulk flow velocities are a few hundreds 
of km/s at scales r < 300 Mpc h~^ (Figure 1), 
the fllters in the WMAP channels do not have 
the sensitivity to measure bulk flow velocities at 
these scales. The differences between the WMAP 
and the Planck filters is expected: the error due 
to instrument noise is largely reduced because the 
noise levels of Planck are much lower than WMAP; 
the wider frequency coverage in the SZ sensitive 
regime also allows filters in the Planck channels to 
suppress the thermal SZ bias. 

10. Conclusion and Discussion 

We investigated Planck performances in de- 
termining bulk flow velocities using the kinetic 
Sunyaev-Zeldovich effect. We characterize the sen- 
sitivity of the bulk flow velocity measurement us- 
ing simulated Planck data combined with three 
representative all sky galaxy cluster surveys: the 
archived ROSAT All-Sky Survey, the Planck clus- 
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RASS 

_'^<^,-^ noise and CMB + thermal SZ 1=30, b= 280 + 

1= 0, b= (X) X 

1= 0, b= 90 ( y ) <> 
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Bulk flow velocity {i<m/s) 
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1=0, b=o(x) X 
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Bull< flow velocity (km/s) 



10000 





Fig. 9. — The 95% upper limit to the error in the angle measurements for all the three cluster surveys. The 
I = 30°, b = 280° data points correspond to the direction of the input bulk flow that the simulations in 
Figure 7 are based on. 
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ter catalog and the cRosita All-Sky Survey. We 
employ two different types of filters, a matched fil- 
ter (MF) and an unbiased kinetic SZ filter (UF), 
to maximize the cluster signal to noise ratio. 
The main results are (see also Table 3): 

1. The use of simulated Planck sky maps in- 
stead of WMAP ones in combination with 
the RASS catalog reduces the velocity that 
can be detected at 95% CL, from ^ 5000 
km/s (~ 10,000 km/s) to ~ 1000 km/s 
(w 1000 km/s) when filtered by the MF 
(UF). 

2. Using all clusters with z < 0.5 we find that a 
bulk flow of 500 km/s would be measured in 
a strongly biased way with the RASS sample 
with both filters. The same velocity is mea- 
sured with 0-5% bias with EASS and Planck 
clusters, with uncertainties of 25-75 km/s. 
These numbers should be compared with the 
rms velocity out to 2; = 0.5 in ACDM: 30 
km/s (95% upper limit: 48 km/s). 

3. If the bulk flow is consistent with ACDM 
prediction of t; = 30 km/s at z = 0.5, our 
analysis pipeline would obtain a 95% upper 
limit to the recovered velocity (when the UF 
is used) oi V = 470 km/s for RASS clus- 
ter, V = 160 km/s for Planck cluster sample 
and t; = 60 km/s for EASS cluster sample. 
This allows us to measure the departure from 
ACDM if the measured bulk flow is in excess 
of these estimates. 

4. Planck can also constrain a recovered bulk 
flow direction. For Vlnput = 500 km/s, the 
95% upper limit to the angle errors (when 
the UF is used) are Aags « 30° for the RASS 
clusters, Aags « 15° for the Planck clus- 
ters, and Aags « 5° for the EASS clusters. 
These uncertainties are lower than the dis- 
crepancies observed in bulk flow directions at 
different optical depth up to now (Figure 1). 
The errors in the directions contained within 
the galactic plane are larger than in the per- 
pendicular one due to the galactic cut at low 
latitudes. 

5. The error in the kinetic SZ dipole is dom- 
inated by the unfiltered CMB and instru- 
ment noise and intrinsic scatter of the Y- 
parameter, with a systematic thermal SZ 
bias. The contamination from extragalactic 
radio sources is negligible and can be safely 
ignored. The noise level is lowest for the 



EASS sample since it contains the most clus- 
ters. 

6. The UF is more sensitive and effective than 
the MF in suppressing the bias induced by 
the thermal SZ in the velocity reconstruc- 
tion. This is in contrast to the performance 
of the two filters on the WMAP maps, where 
the MF was much more sensitive than the 
UF. Differences in performances are to be 
expected, as the two experiments have dif- 
ferent characteristics. As Planck has much 
lower noise, the tSZ contamination now plays 
a more major role in setting the bias and 
errors. Planck increased frequency coverage 
and resolution, however, enables the filters 
(UF in particular) to mitigate the effects 
of the thermal SZ signal, especially when a 
large cluster sample is considered. 

We have not considered here the effect of in- 
frared point sources, correlations in the signal be- 
tween CMB and clusters, large-scale (non bulk 
flow) peculiar velocities, and error determination 
in extracting optical depth from the survey in 
hand. Some of these could present further im- 
provements to this initial study, which demon- 
strate the superior potentials Planck has in de- 
termining the bulk flow. 
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Table 3 
Summary of the results 





noise components 


RASS 


Planck 


EASS 






MF 


UF 


MF 


UF 


MF 


UF 


^500 


noise + CMB + tliermal SZ bias 


1085 ± 82 


720 ± 94 


530 ± 126 


516 ± 73 


507 ± 33 


500 ± 24 


{vx) (75 km/s) 




129 ± 94 


131 ± 138 


108 ± 133 


88 ±68 


80 ±40 


75 ± 29 


{Vy) (-426 km/s) 




-75 ± 105 


-434 ± 133 


-432 ± 123 


-431 ± 73 


-427 ± 37 


-428 ± 26 






1065 ± 82 


526 ± 66 


240 ± 91 


255 ± 49 


256 ± 28 


243 ± 19 


A „500 




62° 


34° 


34° 


14° 


9° 


4° 


iTlos V.binl 




0.71 


0.40 


0.12 


0.09 


0.09 


0.09 


V",bin2 




0.28 


0.14 


0.09 


0.04 


0.02 


0.01 


0"]^Qg \/ bin3 




0.08 


0.05 


0.04 


0.02 


0.01 


0.01 


Clog V,bin4 




0.02 


0.03 


0.03 


0.02 


0.01 


0.01 


Clos V, whole 




0.42 


0.24 


0.08 


0.06 


0.05 


0.05 


Ktr((Koc)) 

(vx) (75 km/s) 


only thermal SZ bias 


1100 ± 60 


724 ± 82 


530 ± 115 


515 ± 63 


511 ± 30 


504 ± 14 




125 ± 73 


142 ± 98 


82 ± 109 


82 ± 60 


104 ± 30 


101 ± 15 


(uy) (-426 km/s) 




-71 ± 75 


-445 ± 110 


-437 ± 112 


-436 ± 63 


-427 ± 34 


-428 ± 14 


(u^) (250 km/s) 




1082 ± 63 


534 ± 46 


253 ± 86 


250 ± 46 


255 ± 29 


244 ± 11 






62° 


33° 


32° 


14° 


9° 


3° 


"■log V.binl 




0.68 


0.40 


0.12 


0.10 


0.07 


0.04 


'^log V,bin2 




0.27 


0.15 


0.06 


0.03 


0.03 


0.02 


flog V,bin3 




0.09 


0.06 


0.04 


0.02 


0.02 


0.01 


flog V,bin4 




0.03 


0.03 


0.03 


0.02 


0.01 


0.01 


flog V, whole 




0.41 


0.24 


0.08 


0.06 


0.05 


0.03 



Note. — This tabic lists the main results of the bulk flow measurements for clusters within z < 0.5 in the three cluster surveys considered 
and for maps containing different components of uncertainties. The first five quantities are based on 100 realizations with input velocity of 
500 km/s at / — 280° and b — 30°. Vj-cc*^' ^^^^ average value of the recovered velocity, (v^) , (v^) , {v^) arc the average of the individual 
directions. For ideal recovery, they should have values (75, -426, 250) km/s. Aqq^^ is the 95% upper limit to the angle error. The last five 
quantities flog v,biiii are the deviation parameters as defined in equation 18, where bin i refers to the binned velocity range considered in the 
calculation: bin 1= 100-500 km/s, bin 2= 500-1000 km/s, bin 3= 1000-5000 km/s, bin 4= 5000-10000 km/s. The whole range is for velocities 
100-10000 km/s. 
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A. Derivation of Filter Kernel 



A.l. Data Model 

In real space, the observation field can be described as 

s^e) = y{0){U - V)B,{e) + n,{e) (Al) 

where j/(0) is the y parameter at the cluster center (defined in equation BIO), f„ = x{e^ + ^)/{e^ — 1) — 4, 
and X = (hu) / {ksTcMB) is the frequency response of the tSZ signal, V = VrmgC/kBTe gives the signal in 
kSZ effect, B^{9) is the convolved cluster profile given by: 



B,{e) = j dQ.'p{e')K{e - e') 

oo 

= ^Bio,.r,O(cos(0)) 

1=0 



where S,o,^ = ^J (An) / {21 + l)biQ, vPio,v^ piP) and 61/ (0) are the cluster spatial profile and beam function 
respectively, and n,y is the noise consisting of instrumental noise and CMB . The noise map and the cross 
power spectrum C of each of the two components satisfy 



A.2. Filtered Field 

Let and be the filter and filtered field respectively at frequency v, then 



Im 



(A2) 



with Uirn,i, = y 2lfT^''"''^*'0^'' 

The total signal in all frequency channels is u{9) = J2u'^'j{(^) 
A useful quantity to consider is the variance of the filtered map: 

al = {{u{e) - {u{e))r) 

We introduce vector notation for the set of frequency dependent quantities, e.g. Bi = {Bi^^^ , Bi^i,^ ,-Bi^vnj} 
Since {u{9)) = udFi — VBi) + ni ^n, where we have used the notation Fi = {fvBio^v}^,, then we 



have 



iVcFi - ycVBi)Si + ni^i - {yA - ycVBi)Si 

. I 



Therefore al = J2i *f Ci$;. 
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A.3. Unbiased Matched Filter for kSZ signal 

In order to derive the optimal filter for the kSZ signal, we want to minimize cr^ subject the the following 
constraints: 

1. The filter is an unbiased estimator of KSZ signal at the source location, such that 

E^'*T = 1 (A3) 



2. The filter should remove the TSZ signal at the source location, i.e. 

J2ml = o 



(A4) 



The functional variation of with respect to subject to constraints A3 and A4 can be obtained using 
Lagrangian multipliers. The Lagrange function is defined as: 



L = J2 *l + Ai(l - E Bi<^f) 



Minimizing the Lagrange function with respect to the filter function 



^ CiSi = XiBi - X2F1 



The job now is to find the constants Ai and A2. Using constraint A3 and A5, 

Y.Bi^I'^i = E [^I'^iBfC^^Bi - X2^iBfCi^Fi 

I I 

^ 1 = Ai E BTC^'B, - A2 E BfC^^Fi 

=> 1 = Ai7 - X2/3 

Similarly, using constraint A4 and A5, 

Y^Pi^I^i = [Ai^/if 'Cf - X2^iFfc;^Fi 
I I 

^ 1 = Ai 5^ Ffc^^Bi - A2 E FfCy'Fi 
^ 1 = Ai/3 - A2Q: 
where a = Cf ^F,, /3 = Ffc^^Bi, and 7 = J], BfQ-^B, 



Solving A6 and A7, we obtain Ai 



OiJ—p 
$UF 



It is easy to verify that this filter kernel satisfies the two constraints. 
The filter can be interpretated in this way: 



.1 

and A2 = ^^"^2 . Thus we get the filter kernal: 



:{aBi-m) 



(A5) 



(A6) 



(A7) 



(A8) 
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1. A = a7 ~ is the normahzation such that ^ Bf^i = 1 

2. The term ^S; comes from the kSZ constraint, ensuring that the filter gives the kSZ signal at the source 
location. 

3. The term —-^Fi comes from the TSZ constraint. Its purpose is to suppress the tSZ signal such that 

tSZ signal vanishes at the source location. 

4. The filtered fields at each channel are weighted by the inverse of the covariance matrix Cf"*^ and then 
combined to form the final filtered signal. 

A.4. Matched Filter 

The derivation of the matched filter is similar to that for the unbiased matched filter, except that con- 
straint A3 is used. Therefore we have the Lagrange function: 

I I 
^ = E [Ci^i - ^Bi] = 

^ = C^^iXBi) (A9) 
Using constraint A3 and solving for A, we find A = I/7. Thus we get the filter kernel: 

$MF ^ -J—Bi (AlO) 
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B. Simulations of full-sky SZ maps 

B.l. Mass and Redshift Distribution of the cluster sample 

The number of galaxy clusters per mass and redshift bin can be estimated using the comoving mass 
function: 

dN . ^ . . ^ dn . dV , . 
, {M,z) = An-—(M,z)——(z) 

dMdz^ ' ' dM^ ' 'dzdn^ ' 

where AfJ is the soHd angle in a given direction, and -^{M,z) is the mass function which describes the 
number density of galaxy clusters per mass bin at a given redshift z. Here we adopt the Jenkins mass 
function (Jenkins et al. 2001) with the fitting formula given by: 



= 0.315 exp(- |lnCT-i +0.64|^-^*) 

where pm is the present matter density, a is the variance of the mass fluctuation cr^ = k^dkPm {k) \ Wu (fc) | , 

where Pm{k) is the matter power spectrum P{k,z) = Ak'^T'^{k,z)D'^{z), 

Wc use the transfer function T{k, z) from Eiscnstcin & Hu (1998) that accounts for all baryonic effects in 
the matter transfer function on the large scale and an improved version of the BBKS transfer fitting formula. 
The growth function is defined as D+{z) = Di{z)/Di{0). We take the approximated form from Lahav et al. 
(1991) and Carroll et al. 1992: 

Di{z) = (1 + z)-i^^Mf) [0„z4/7 _ ^^(^) ^ (1 ^ n^z/2){l + r!A(z)/70)] (B2) 

By integrating over the mass and redshift interval of interest, we obtain the number of clusters N in a 
redshift and mass bin: 

NamAz=^tt f dz [ dM^dz{M,z) (B3) 
Jaz Jam "J"' 

We then assign TV' clusters to a bin (AM, Az) from a random Poisson distribution with an average 
value given by equation B3. The mass range is taken to be 10^'* — IO^^Mq with a bin size of log(AAf) = 
0.1. The minimum mass of cluster at redshift z that can enter the simulated catalog is given by Mmin = 
max[10^^MQ, Miiin], where Mii,„ is the limiting mass of the corresponding survey from equation 9 and 10. 
The redshift range z is — 0.3 with a bin size of Az = 0.02. The number of clusters from this halo model is 
a strong function of as shown in Figure 4. 

Having obtained the mass and redshift distributions of the clusters, we distribute the clusters over the 
whole sky. For simplicity, we ignore the spatial correlation among clusters and assign a random position to 
each cluster. 

B.2. Modeling SZ signal of individual clusters 

The SZ signal of galaxy clusters is characterized by the Compton y parameter, which in turn depends on 
the gas properties of individual clusters. The next step is to model the electron density ne(r), temperature 

Te{r) and the Integrated Compton Y parameter. 

B.2.1. Gas properties from scaling relations 

Assuming that clusters are self-similar, virialized and isothermal, we obtain the mass-temperature relation 
from the virial theorem (Kaiser 1986): 



1/3/ ^ x2/3 



keV (B4) 
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where /3t = 0.75 is the normahzation constant under the assumption of hydrostatic cquihbrium and isother- 
mality, Ay is the mean overdensity of a viriahsed sphere which we have calculated using the fitting formula 
from Pierpaoli et al. (2001), Md is the mass of the cluster within the virial radius r^ir, = An/S AvPcrit?"vir- 
Solving for Tyir gives, 



9.5103 /a„Av(z)\ '''V Afci ^ 



1/3 



(B5) 



The assumption that clusters arc isothermal is in good agreement with XMM-Newton observations of outer 
cluster regions which find that the cluster temperature profiles are isothermal within ±10% up to » rvir/2 (Ar- 
naud et al. 2004). 

The electron density within the virial radius described by a ^ model ( Waizmann & Bartelmann 2009 and 
reference herein), 

/ ^2 \ -3/2/3 

ne(r) = rieo I 1 + ^ 



(B6) 

where Vc is the core radius. Geisbiisch et al. (2005) obtain a relationship between the core radius and the 
virial radius: 

r,(z) = 0.14(1 + z)i/Vvir (B7) 

We choose /3 = 2/3 to match Waizmann & Bartelmann (2009) which describes the X-ray surface brightness 
profile of observed clusters. 

The central electron density is normalized by JJ^""" ne{r)dV = N^, in which A^e is the total number of 
electrons within the cluster given by: 



1 + ./H 

2mD 



/gasMcl 



(B8) 



where /h is the hydrogen fraction, /gas = ^h/^m is the baryonic gas mass fraction of the total cluster mass. 
Here we take /gas = 0.168 from WMAP 5 year results and /h = 0.76. We find. 



A^. 



?^eO = 



C 1(0.14(1+2)1/5) 



■tan-i(0.14(l + z)i/5) _ | 



(B9) 



The Comptonization parameter is proportional to the product of the electron temperature and density, 
and is given by y = / dl{kBTe) / {meC^)neCrT ■ Performing this integral along the line-of-sight of the cluster, 
we obtain: 



(0.14(l+z)^ 



5? 



1 + ^ 



(BIO) 



where 9c = rd Dp^. The SZ signal is then calculated using this compton y parameter for each cluster within 
the virial radius. It is easy to show from equation BIO that y drops to zero when 9 > ^vir- The analytic total 
Comptonization parameter Y is obtained simply by integration equation BIO over the cluster size: 



Y = [ d^y{9) = D^^z)^^ [ dVn,T, = 



where Da is the angular diameter distance to the cluster. We then scale the pixels within the disc of angular 
size ^vir by a factor N = ^ y^npi > where is the y value in the center of pixel i and rfOpix is the pixel solid 
angle. 

There is an expected cluster mass-optical depth scatter in nature due to the fact that unknown cluster 
physics can affect the exact form and normalization of the scaling relations, e.g. Y — M relation. Sehgal et al. 

(2007) have estimated the typical intrinsic scatter in the Y — M relation with simulations of the ACT clusters 
and found that this is about 15% in Y. While ACT clusters represented in those simulations have masses 
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that are lower than those expected in our samples (in particularly Planck), we do not expect bigger objects 
which are more virialized to have a larger scatter. We therefore introduce a 15% scatter in the Y parameter 
in our simulations for all cluster samples. In principle, there is also a measurement error for optical depth 
due to inaccuracy of SZ photometry. However, we ignore this effect because such uncertainty is not yet well 
characterized for Planck. 
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